Nuclear matter properties in relativistic mean field model with ct-uj coupling 
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The possibility of extending the linear a-uj model by introducing a a-u) coupling phenomenolog- 
ically is explored. It is shown that, in contrast to the usual Walecka model, not only the effective 
nucleon mass M* but also the effective a meson mass m* and the effective oj meson mass are 
nucleon density dependent. When the model parameters are fitted to the nuclear saturation point 
(the nuclear radius constant ro = 1.14fm and volume energy a\ — 16.0MeV) as well as to the effec- 
tive nucleon mass M* = 0.85M, the model yields m* — 1.09mo- and m* — 0.90mui at the saturation 
^— ^ ', point, and the nuclear incompressibility Kq — 501MeV. The lowest value of Kq given by this model 

■ by adjusting the model parameters is around 227MeV. 

o : 
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. As the starting point for the microscopic relativistic description of nuclear many-body system, within the framework 
' of quantum hadrodynamics, the well-studied linear a-cu model has been proved to be able to describe the saturation 
t . and other properties of nuclear matter However, this model yields the nuclear incompressibility Kq around 550MeV 
which is unacceptably high, and also the effective nucleon mass M* around 0.54M which seems uncomfortably low. 
In order to remedy this situation, many works have been done on the extension or generalization of this model. 
Further inclusion of nonlinear self-interaction of a mesons has been proved to be successful in reducing the nuclear 
incompressibility Kq significantly, in spite of some controversy about the stability of the system specified by the fitted 
parameters On the other hand, the density-dependent N-cr interaction given by the Zimanyi-Moszkowski model 

^— I 1^ is shown to be able to give reasonable nuclear incompressibility Kq and effective nucleon mass M* as well. Since 
, the Zimanyi-Moszkowski model contains a coupling between a and uj mesons, it is interesting to investigate what is the 
' direct effect of the coupling between a and lu mesons, when detached from the density-dependence of the N-cr coupling. 
The purpose of this letter is to explore the possibility of extending the linear a-uj model by introducing a a-uj coupling 

■ phcnomcnologically, and to investigate whether this coupling alone is able to reduce the nuclear incompressibility Kq 
"*7^ [ and in the same time to increase the effective nucleon mass M*. 

'"^ ■ Let us start with the following Lagrangian density: 

§ : C= ^[j^iid^ - g^ion - (M - 3,0)]^ 

: 



> 

in 



-^Ff^F^ + ^mlio^u:', , = !imdl, (1) 

where F'^" = d'^uj'^ — d'^uj'^, -0, and cu are the nucleon, a and u meson fields with masses M, rria- and rriuj, respectively, 
while g^, g^j are the respective coupling constants, and g'^ is the <t-uj coupling constant, rj — i gives a a-uj coupling 
similar to that introduced by usual covariant derivatives, while rj — 1 gives a gauge- like virtual coupling similar to 
the imaginary coupling [Q. Mathematically, it is equivalent to add a term rj^g'^ui^j^ui^^ to m^, which is the same as the 
rj2 term introduced in the nuclear effective field theory |^ . The present model reduces to the original Walecka model 
when the a-uj coupling constant is zero, = 0. 

For static nuclear matter, the field equations derived from this Lagrangian are reduced to the following equations 
in the mean field approximation: 

{ird^-g^fLOo~M*)^ = Q, (2) 

Tn*^'P = (3) 
■m*J'ujQ = gujRN, (4) 
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where the effective nucleon mass M*, effective a meson mass m* and effective w meson mass m* are defined respec- 
tively as 

M* = M - g^d^, (5) 

*2 2,2 /2 2 /n\ 

*2 2 2 12 i2 f'7\ 

and ps — {t/ji/j) is the scalar density, pn = (V'7*''0) the baryon density. For calculating these effective masses, the 
above three equations can be rewritten as the following self-consistent equations: 

> M* a Ps 



TO 



TO^ 

where po is the standard nucleon number density, a, Ua and are the dimensionless composite parameters defined 
respectively as 

m^M m^jm^ rrv^m^ 

It should be noted that and are positive for rj = \ while negative for = i. 

The nuclear matter equation of state derived from Lagrangian density (|l]) can be expressed in terms of the nuclear 
energy density £ as e — £ / pN — M , and 

£^£k+£a+£u., (12) 

£k = ^y2 F,{k,/m, (13) 

i—p,n 

£. = lil-OMps, (14) 

£. = lii^OyMp., y-^-^^^^^, (15) 

2 s(l — u) 

where kp and fc„ are the proton and neutron Fermi momenta respectively, and the function Fm{x) is defined as 

F„,{x)= ( dxx'^"'^/l + x^. (16) 

^0 

The baryon density pN and scalar density ps can be expressed as 



Ps 



J2 flih/m, (18) 
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where the function fm{x) is defined as 



dx- 



10 VT+x^ 

Having the equation of state, the pressure p can be derived as 

OpN o S 



(19) 



(20) 



Instead of Fermi momenta kp and fc„, we will use the nucleon density pjy = Pn + Pp a-nd relative neutron excess 
S = {pn — Pp)/ Pn as independent variables for describing the nuclear matter Q. At the standard state pjv = Poj 
(5 = 0, the pressure should be zero, 



P(po,0) = 0. 

In addition, the depth of the equation of state is related to the nuclear volume energy ai as 

e(po,0) = -fli. 



(21) 



(22) 



The standard nucleon number density po defined by Eq.(pTl) is related to nuclear radius constant tq and Fermi 
momentum kp oi standard nuclear matter as po = 3/47rrQ = 2fc|,/37r^. 

The generalized nuclear incompressibility K{pN,S) can be defined Q and derived as 



K = 9 



dp 3{p + £) + U£^ 
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1 — 2v V dpN 
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2{l^v) dS, 
^-^l^vP^^d^ 
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(25) 



dps ^ , ^.Pn 
PN^r~ = Q + 3(Ps - Q)- 
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(27) 



i—p.n 

The usual nuclear incompressibility Kq can be obtained from K(pN^ S) as 

ro = i^(po,0) = 9(p^|^)^, 



(28) 



(29) 



the subscript stands for the standard state pN = po and 5 = 0. In addition, the following formula for symmetry 
energy J can be derived: 



~ 2 9^2 



1 



G^kiTMm 



(30) 



Around the standard state pN = Po, S = 0, the nuclear equation of state and thus the nuclear matter properties are 
specified essentially by the standard density po, volume energy ai, symmetry energy J, incompressibility Kq, density 
symmetry L and symmetry incompressibility Kg p], where L and Kg are, respectively. 
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For the description of nuclear matter, there are three independent parameters a, and a^j in the present model. 
The saturation point of nuclear matter (po,ai) defined by Eqs.(|T]) and ( p^ can be used as input to determine 
and auj as function of a. The way is as follows. At the standard state (pjv = po,S = 0), these two equations can be 
solved for s and y, 

3(1 -OMp. (33) 



3po{M -ai)-2£k- ^Mp, 

ii) — 2£k — 
(1 - OMP. 



2po{M ~ a,) - 2£k - (1 - QMp, 

n ?V777 ' (34) 



and thus from Eq.(|15| 



sy + s — 1 
sy + 2s-2' 



(35) 



For given po and ai, s and v can be calculated by Eqs.(|33|)-(|35|) for a chosen ^. Then a, and a^, can be calculated 
by Eqs.(||)-(^o|) for a chosen In our calculation, the nuclear radius constant tq = 1.14fm which corresponds to 
the standard nucleon number density po = O.lGlfm^'^, the nuclear volume energy ai = 16.0MeV, nucleon mass 
M = 938.9MeV and constant Tic = 197.327053MeV-fm are used. 

Fig.l plots a„ and auj as function of a. At the point a = aw = 0.4908 (indicated by an arrow), we have 
OLa = OL^ — 0, and the present model reduces to Walecka model. For a > aw, 0£a and a^j are negative, which 
corresponds to rj = i. For a < aw, ol^ and a^ are positive, which corresponds to r\ — 1. There is a lower limit 
curain = 0.1037, bclow this poiut there is no physical solution, as the effective lo meson mass becomes imaginary, 
^5 < 0. Fig. 2 shows the dimcnsionless effective nucleon mass ^, the dimcnsionless effective a meson mass and the 
dimcnsionless effective uj meson mass as function of a. At the point aw , £, = 0.5437 and ^o- = = 1. From aw to 
the right, ^ decreases as a increases, and we have <^ct < 1 and > 1. From aw to the left, ^ increases as a decreases, 
and we have ^cr > 1 and < 1. Fig. 3 gives Kq, J, L, and Kg as function of a. At aw, Kq — 553MeV, J = 20.2MeV, 
L = 70.6MeV, and Kg = SSMeV. As a decreases, Kq decreases at first and then increases slowly, passing through a 
very flat plateau, and finally decreases rapidly to a lower limit 227MeV. J and L as well as Kg decrease slowly as a 
decreases. It is worthwhile to note that Kg is negative in the low a side, in opposition to what is obtained in the usual 
a-Lu model. Experimentally, Kg obtained from the isoscalar giant monopole resonance energy is between — 566± 1350 
to 34 ± 159MeV On the other hand, Kq increases to very high values rapidly as a increases. In this aspect, the 
case of 77 = i is not acceptable. 

The dimcnsionless effective nucleon mass ^ at the saturation point can be used as the third input to fix the parameter 
a. For example, if ^ = 0.85 at the saturation point, we have a — 0.1822. Correspondingly we have Ofo- = 0.1134, 
a^ = 0.2938, = 1-085, = 0.895, Kq = 501MeV, J = 13.8MeV, L = 30.1MeV, and Kg = -38MeV, where the 
values of and are those at the saturation point. Instead of ^, the nuclear incompressibility Kq could be also 
used to fix a. If Ka = 250MeV is chosen, it yields a = 0.1050, a„ = 0.00783, a^ = 2.642, ^ = 0.939, = 1-295, 

— 0.328, and Kg = — 31.2MeV, where the values of ^, ^o- and are those at the saturation point. 

The symmetry energy J can be increased by including the p meson field in the present model. In this case, the 
additional term in the Lagrangian density is 

Cp = -^B^^.B^"^ + irn^b^.b^ - g^^j^T-h^^^P, (36) 

where B'"' = d'^h'^ — d'^h'^, is the p meson field with mass nip and coupling constant gp, and t are the isospin 
matrices. The resulting nucleon field equation, in the mean field approximation for static nuclear matter, has an 
additional term —gpT^^^h^: 

ill'' dp - g^fiOQ ~~ gpT3-fX - M*)4, = 0. (37) 
The equation for the p meson field is 
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&o = ~^, (38) 

while Eqs.(|) and (|) are the same. The additional terms to the nuclear energy density, pressure, generalized incom- 
pressibility, symmetry and density symmetry energies are respectively 

£p = l-apMp%6'^/po, Pp = £p, Kp = ^ -^^pM, Lp = 3Jp, (39) 



where 



mjM 



(40) 



is the p meson dimensionless composite parameter. It is easy to prove that there is no contribution to Kg from the p 
meson field. 

It can be seen that, as £p = at the standard state where 6 = 0, the inclusion of p meson produces no change 
in nuclear energy density £, pressure p and generalized incompressibility K at the standard state. Therefore, the 
parameters a, and a^j have no change, and thus the dimensionless effective masses ^, ^o- and are the same even 
the p meson field is included. The parameter ap can be fixed by using the symmetry energy J as the fourth input, if 
parameters a, and a^^ have been fixed. 

In conclusion, the possibility of extending the linear a-uj model by introducing a a-u coupling phenomenologically 
is explored. It is shown that, in contrast to the usual Walecka model, not only the effective nucleon mass AI* but 
also the effective a meson mass m* and the effective co meson mass m* are nucleon density dependent. When the 
model parameters are fitted to the nuclear radius constant rg = 1.14fm and volume energy ai = 16.0MeV as well as 
to the effective nucleon mass M* = 0.85M, the model yields m* = I.OOtoo- and m* = 0.90m^ at the same nuclear 
saturation point, and the nuclear incompressibility Kq = 501MeV. This incompressibility seems too high. On the 
other hand, if the model parameters are fitted to Kq = 250MeV, in addition to ro = 1.14fm and ai = 16.0MeV, it 
yields M* = 0.938M, m* = 1.295mcr and m* = 0.328m(j at the nuclear saturation point. This effective lu meson mass 
seems too low. In addition, the lower limit of incompressibility Kq = 227MeV is not low enough. Therefore, even this 
model is able to reduce the nuclear incompressibility and also to increase the effective nucleon mass simultaneously, 
there is still not enough degrees of freedom to adjust the parameters in order to give more reasonable result. In this 
aspect, additional physics should be included for improving this model. 
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FIG. 1. Qct and a^, as function of a. The point aw = 0.4908 indicated by an arrow corresponds to the usual Walecka model 
where acr ~ a^i = 0. 



FIG. 2. The dimensionless effective nucleon mass ^, the dimensionless effective a meson mass ^o- and the dimensionless 
effective ui meson mass as function of a. At the point aw = 0.4908 indicated by an arrow, ^ = 0.5437 and = ^uj ~ 1- 



FIG. 3. Ko, J, L, and Ks as function of a. At the point aw = 0.4908 indicated by an arrow, Ko = 553MeV, J = 20.2MeV, 
L = 70.6MeV, and Ks = 88MeV. At the lower limit Q,„i„ = 0.1037, Kq = 227MeV. 
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